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Abstract 

We study the line gamma ray signature from the Galactic Center originating from Bino-like dark 
matter pair annihilation, when CP is not conserved. We consider the case that the Bino is nearly 
degenerate with a sfermion in mass. Although an indirect detection of Bino-like dark matter is 
difficult in general, the signal of Bino-like dark matter pair annihialtion in this case is enormously 
enhanced due to the threshold singularity. 
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Recent observations of cosmological and astrophysical quantities determine the mean 
density of matter and baryons in the universe precisely, and the existence of non-baryonic 
dark matter has now been established [lj]. However the composition of the dark matter is 
still an unresolved problem. Weakly interacting massive particles (WIMPs) are considered 
to be a prominent candidate for dark matter Q]. 

The minimal supersymmetric standard model (MSSM) provides the stable lightest su- 
persymmetric particle (LSP), if R-parity is conserved. The neutralino LSP can be the dark 
matter due to the nature of the WIMpj^]. The neutralino is a linear combination of the 
supersymmetric partners of the U(l) and SU(2) gauge bosons (Bino and Wino), and those 
of the Higgs bosons (Higgsinos). 

In many breaking scenarios, the Bino-like neutralino (B) is predicted as the LSP. The 
viable models with bino-like dark matter require the presence of mechanisms which suppress 
the density. One of the mechanism is the coannihilation|4]. It works when there are other 
particles whose masses are nearly degenerate with that of the dark matter. In this case, 
the relic abundance of the dark matter is determined not only by its own annihilation but 
also by that of the slightly heavier particles, which will later decay into the dark matter. 
For the Bino-like dark matter, the candidate for the coannihilating particle is the sfermion 
/, which is the supersymmetric partner of the fermion /. One such particleis the stau f, 
and the coannihilation between Bino and stau has frequently been discussed in the past 
works, especially in the minimal supergravity model (MSUGRA). Another example is the 
stop t. This possibility is considered in the context of electroweak baryogenesis in MSSM. 
The existence of the Bino-like dark matter having (9(100) GeV mass is predicted in both 
cases. 

Futhermore, the effects of CP violation have been considered not only on the flavor 
physics but also on the dark matter physics. In dark matter physics, the sfermion sector is 
extensively studied 0. 

For investigating the nature of the dark matter, many ideas for its detection are proposed 
and some of the experiments are now operating. Among those, the detections of exotic 
cosmic ray fluxes 6], such as gamma rays, positrons, anti-protons and neutrinos from dark 
matter annihilation are feasible techniques. In particular, an excess of monochromatic (line) 
gamma rays originating from the dark matter pair annihilation to two gammas would be a 
robust evidence if observed, because diffused gamma ray background must have a continuous 
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energy spectrum. The satellite detectors and the large Atmospheric Cerenkov Telescope 
(ACT) arrays may search for the exotic gamma rays from the Galactic Center Q]. 

In this letter, we study the line gamma ray flux from the annihilation of the Bino-like 
dark matter which is nearly degenerate with stau or stop in mass and has CP violating 
interaction with those particles. Interaction of the Bino-like neutralino is so weak that the 
line gamma ray flux from the annihilation has been considered to be very small without 
CP violation. We find, however, that the annihilation cross section is enormously enhanced 
due to the threshold singularity when CP is violated compared to that of the CP conserved 
case. Since the dark matter is non-relativistic (the typical velocity is O(10~ 3 )), the relevant 
state for the incident Bino pair is 1 S'o, which is CP odd due to the Majorana nature, while 
the sfermion pair in the 1 S'o state is CP even due to the bosonic nature. When CP is not 
conserved, the transition from the incident Bino pair to a sfermion pair in the 1 S'o state 
takes place without the velocity suppression v in the non-relativistic limit, and the sfermion 
pair can annihilate into two gammas. 

The source of the CP violation is contained in the off-diagonal element of the mass matrix 
for the third generation sfermion M~, 



Ml + m 2 f + m 2 z cos 2/3 (T 3L - Q sin 2 9 W ) 
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1 mJ LR Ml + m 2 f + m| cos 2(3 Q sin 2 9 W 

where 

m /Li? = — m r(A* + //tan/3) , (for stau) , 

= -m t (A* t + li cot (3) , (for stop) , (2) 

Ml and Mr are the left- and right-handed sfermion masses, and tan (3 is the ratio of the 
vacuum expectation values of two Higgs fields. The charge and isospin of the sfermion are 
denoted by Q and T 3L , respectively, and 9w is the Weinberg angle. The parameters rrif, m z 
and \x are the fermion, Z boson and supersymmetric Higgsino mass, respectively. The mass 
matrix of stau or stop is least studied experimentally, in particular the trilinear A T or A t 
parameter is not constrained from stringent EDM (Electric Dipole Moment) experiments. 
Thus, there can be a wide space in A T or A t to accommodate large CP violation. 

For following discussions, we define the sfermion mixing matrix for diagonalizing the mass 
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matrix, U^M 2 -U = diag.fm^ ,m 2 - ), as 
/ x h h 

( h )=u(' l ) = ( (') . (3) 

\fn) \h) \-sm6 f e-V cos 9, J \f 2 J 

where we use the notation that fx has smaller mass than / 2 , m; < m^. Hereafter, fi is 
denoted by / and its mass by for simplicity. We assume that the lightest sfermion / is 
almost degenerate with B in mass, which is needed for the coannihilation process. We also 
assume that other SUSY particles are so heavy, in general much heavier than /, that all the 
experimental constraints are satisfied in the presence of the large A T or A t . 

The transition from the incident Bino pair to the sfermion pair is caused by the Yukawa 
interaction B — f — f. Its strength is given by the gauge coupling due to supersymmetry, 

C tr = f*B (a L P L + a R P R ) f + h.c. , (4) 

where, 

a L = -V2g'{Q-T 3L )cos6 f , (5) 

a R = -V2g'Q sin 6 f e* 7 ' , (6) 

g' is the U(l)y gauge coupling constant of the standard model. The effect of CP violation 
is induced by the relative complex phase e nf . 

The threshold singularity appears in the calculation of the annihilation cross section 
due to the degeneracy in mass between a Bino pair and an intermediate sfermion pair and 
their non-relativistic motions. The intermediate sfermion pair forms a quasi stable bound 
state when the binding energy is approximately equal to the mass difference, and hence the 
annihilation cross section is enhanced. In a Feynman diagrammatic picture, higher order 
terms such as a ladder diagram depicted in Fig^ significantly contribute to the cross section, 
because an intermediate sfermion pair is almost on-shell. Thus, we have to resum the ladder 
diagram for a reliable calculation. The non-relativistic (NR) effective theory is useful to 
treat the singularity [8], which is often used in calculations of the threshold production of 
particles, the quarkonium mass spectrum and so on. The Schrodinger equation for the 
non-relativistic Bino pair annihilation is derived from the NR effective Lagrangian as its 
equation of motion. Solving the Schrodinger equation is equivalent to the resummation of 
ladder diagrams in Fig^ 
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FIG. 1: Ladder diagram contributing to the threshold singularity in BB — > 77 process. The 
character g means gluon, which is important when / is top. 

The derivation of the Lagrangian here is essentially the same as that in Ref . |9| . Therefore 
we show only the result. The NR effective Lagrangian for two-body states, 4>§{r) (— BB/2) 
and 4>f{f) (— /*/), is given as 

C = *t(r) 

where $>*(r) = (0°~(r), ^~(f)), r is the relative coordinate (r = \ f\), E is the internal energy 
of the two-body stated = mv 2 /4, m is the Bino mass, v is the relative velocity of the incident 
dark matter pair (v = \v\). The parameter C is the color factor, which takes a value C = 1 
for stau and C = 3 for stop. The potential V(r) has 2x2 matrix form. The (1,1) component 
describes the force acting between the sfermion pair and the dominant contribution is the 
coulomb force for the stau case, while it comes from the QCD interaction for the stop case, 

Vn(r) = 25m- - , (8) 
r 

where 8m is the mass difference between the sfermion and Bino, d — a for the stau case, 
and d = 4a s /3 for the stop case. We set the fine structure constant a = 1/128, and the 
strong coupling a s is determined by solving the one-loop renormalization group equation at 
the energy Ama s /3; 

a s (4mx/3) = x (9) 

where a s (/x) is the effective coupling at scale \i. Eq. (jHJ) is derived from the consistency 
equation that the typical momentum of the intermidiate gluon in the t-channel is about md, 
when the sfermion pair forms a bound state by the attractive force of the second term of 
right hand side in (JHJ). The (1,2) (= (2,1)) component represents the transition between 



E + —)l- V(f) + iCT5 3 (r) $(f) 
m 



(7) 
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the sfermion pair and the Bino pair, and is derived from the interaction in eq.(J3J, 

V 12 (r) = f -e- m f r , (10) 

where 

/ = -^^- ^{a R a* L ) oc sin 20/ sin 7/. (11) 

Since the transition takes place due to the CP-violating interaction, the component is pro- 
portional to sin7j as expected. The (2,2) component of the potential V is zero, because 
there are no forces acting between Bino pair in the non-relativistic limit. The imaginary 
(absorptive) part of the effective Lagrangian T is 

= 5 — OnOji ■ (12) 

m 2 

The sfermion pair annihilates into two gammas through this part. 

The amplitude for the Bino dark matter pair annihilation is obtained by solving the 
Schrodinger equation 



V 2 

-— + V r) 

m 



$(f) = E&{r) , (13) 
where the boundary condition is given by 

<&(r — > oo) = I , I , k = —v . (14) 




The annihilation cross section of Bino-like dark matter to two gammas is calculated as 

^=($|CT|$) = ^^|0 / (O)| 2 . (15) 

We perturbatively expand the solution of eq. (|T5|) by the off-diagonal element of the poten- 
tial (|TU|). Under the boundary condition ([Tljl. the solution to 1st order in the perturbation 
is given by 

2f f°° e ~m f y/dm 

*' {0) = -il dv ^ tMA <i)- (16) 

where, 

(W = - r (-2^) e ^ lC,(1 -2^' 2 - 2 ^)- < 17) 



where a = -^^{25ra — ^j-), T(x) is Euler's gamma function, and the function U(a,b,c) is 
the confluent hyper geometric function defined by 

U(a } b } c) = — - / dt e~ ct t a ~ l (1 + t) h - a - 1 . (18) 

Then the annihilation cross section of Bino-like dark matter is given by substituting (JTfiJl 
for (|15j). In Fig. the numerical results for the annihilation cross section to two gammas 
are shown. We assume that the sfermion mixing is maximal and CP violating phase is also 




FIG. 2: Contour plot of the annihilation cross section of Bino-like dark matter to two gammas, (a) is stau 
case, and (b) is stop case. We take the sfermion mixing as 9f — 7r/4, and CP violating phase as 7/ = n/2. 



maximal. Our result is sensitive to the sfermion mixing and CP violating phase. However 
this result is maintained as long as the sfermion mixing is not so small and the CP violating 
phase is 0(1). 

The annihilation cross section is enormously enhanced, if 8m is approximately equal to 
the binding energy of the sfermion pair. The binding energy is determined by 

md 2 „ mv 2 , , 

— = 25m-—, 19 
An 2 4 v ' 

where n is a principal quantum number. When n — 1, 5m/m ~ d 2 /8. Indeed the enhance- 
ment is realized there as shown in Fig. El 

We note that the neutralino annihilation cross section to two Grammas has been calculated 
in the CP conserving easefl, and it tate.y exeeeds a few times HT- cnV for Bino-lifce 
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neutralinos [6(. On the contrary, in our case, the annihilation cross section can be enhanced 
by several orders of magnitude compared to the CP conserving case. 

The enhancement of the annihilation cross section gives significant impact on an indirect 
detection for the dark matter search, because the flux of line gamma rays is proportional to 
the cross section. We discuss the line gamma ray flux from the dark matter annihilation in 
the Galactic Center, which is given by 

F ta = 1.9 x 10-" ( ( <">, , ) J(Afi) Ail cm-V , (20) 

where (av) is the thermally averaged annihilation cross section, and AQ is the angular 
acceptance of the detector, taken to be 10~ 3 for our calculation, which is a typical value for 
ACT detectors. The flux depends on the dark matter profile through the parameter J(Afi) 
in eq. ipnjl . 

7= [ — I - ( p Y 

= J An Ml J Vme of site 8Tkp^ 1 0.3 GeVcm-3j ' ' " ! 

where the integral over / is performed along the line of sight. The parameter J(AQ) is 
estimated in various halo model, and takes values 10 < J 10~ 3 ) < 10 6 . We take J(10~ 3 ) = 

n □ 

1352 for the NFW profile which is a moderate value[2j. 

Our prediction for the line gamma ray flux from the Galactic Center is shown in Fig. 01 
The line gamma ray flux is enourmously enhanced, and will be observed by the ACT detec- 
tors. 

To conclude, we have calculated the line gamma ray flux from the Galactic Center origi- 
nating from the Bino-like dark matter annihilation in the CP violating case. We have showed 
that the line gamma ray flux is enourmously enhanced due to the threshold singularity when 
the degeneracy between Bino and sfermion is sufficiently small. 

HESS and CANGAROO-II have already taken data from the Galactic Center. Some 
parameter region may have been excluded by these ACT detectors. However, the results 



12(|, and hence we do not compare our 



of HESS and CANGAROO-II show the discrepancy 
prediction with these results in detail. Our predicted flux is large and hence the line gamma 
ray will be observed by the ACT detectors in the near future, even if the neutralino is 
Bino-like. 
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FIG. 3: Contour plot of the line gamma ray flux from the Galactic Center originating from Bino-like dark 
matter pair annihialtion. (a) is stau case, and (b) is stop case. We take Aft — 10~ 3 , J(10~ 3 ) = 1352, 
Of = 7r/4, and 7/ = n/2. 

Acknowledgments 

We thank M.M.Nojiri for useful discussion. The work of J.S is supported by the Grant- 
in-Aid for Scientific Research on Priority Area No. 16028202 and 17740131. 



[1] D. N. Spergel et al. [WMAP Collaboration], Astrophys. J. Suppl. 148,175 (2003); C. L. Ben- 
nett et al, Astrophys. J. Suppl. 148,1 (2003). 

[2] For reviews, G. Jungman, M. Kamionkowski and K. Griest, Phys. Rept. 267 195 (1996); 
L. Bergstrom, Rept. Prog. Phys. 63,793 (2000); C. Munoz, Int. J. Mod. Phys. A 19, 3093 
(2004); G. Bertone, D. Hooper and J. Silk, Phys. Rept. 405,279 (2005). 

[3] H. Goldberg, Phys. Rev. Lett. 50, 1419 (1983); J. Ellis, J.S. Hagelin, D.V. Nanopoulos, 
K.A. Olive, and M. Srednicki, Nucl. Phys. B 238, 453 (1984). 

[4] K. Griest, and D. Seckel, Phys. Rev. D 43, 3191 (1991); S. Mizuta, and M. Yamaguchi, Phys. 
Lett. B 298, 120 (1993); S. Profumo and C.E. Yaguna, Phys. Rev. D 70, 095004 (2004). 

[5] T. Falk, K.A. Olive, and M. Srednicki, Phys. Lett. B 354, 99 (1995). 

[6] L. Bergstrom, P. Ullio, and J.H. Buckley, Astropart. Phys. 9, 137 (1998) 



9 



[7] A. Morselli, SLAC-REPRINT-1997-111; K. Tsuchiya et al, Astrophys. J. 606, L115 
(2004)[CANGAROO Collaboration]; J. A. Hinton et al, New Astron. ReV. 48, 331 

(2004) [HESS Collaboration]; C. Baixeras et al, Nucl. Phys. B,Proc. Suppl. 114, 247 
(2003) [MAGIC Collaboration] 

[8] W.E. Caswell and CP. Lepage, Phys. Lett. 167, 437 (1986). 

[9] J. Hisano, S. Matsumoto, M.M. Nojiri, Phys. Rev. D 67, 075014 (2003); Phys. Rev. Lett. 92, 
031303 (2004); J. Hisano, S. Matsumoto, M.M. Nojiri, O. Saito , Phys. Rev. D 71, 063528 

(2005) . 

[10] L. Bergstrom, and P. Ullio, Nucl. Phys. B 504, 27 (1997); Z. Bern, P. Gondolo, and M. Perel- 
stein, Phys. Lett. B 411, 86 (1997); P. Ullio, and L. Bergstrom, Phys, Rev. D 57, 1962 
(1998). 

[11] J.F. Navarro, C.S. Frenk, and S.D.M. White, Astrophys. J. 490, 493 (1997). 
[12] F. Aharonian et al. [HESS Collaboration] , |arXiv:astro-ph/0408145| 



10 



